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For appropriate parameters, the ground state for the Schrodinger and Ampere coupled equations 
in a cylindric domain does not have axial symmetry. 
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Symmetry breaking — a situation in which the solution 
of a problem has lower symmetry than the problem it- 
self — is usually a room for interesting physics (e.g. jij). 
It is therefore gratifying to note that symmetry break- 
ing is found in one of the most common problems we 
encounter: that of a condensate of particles interacting 
with a magnetic field, in a time independent state. 

Our model is as follows. We define a "thermodynamic 
potential" density 



3=|Vx Ai| 



(*V - A; - br9)^\ 2 



(1) 



where ip is a wave function, A; the magnetic potential 
induced by ip, b is the external magnetic field, which is 
taken as uniform, and r, 6 are cylindrical coordinates, 
with the z-axis parallel to the field, b and A; are nor- 
malized so that no coefficients are required in (jl]). The 
sample in which electrical currents can exist will be an 
infinitely long cylinder, and therefore VxAj = outside 
of the sample. The first term in (|lj) gives the contribu- 
tion of the magnetic field to the thermodynamic poten- 
tial, subject to the constraint of a fixed external field; the 
second term may be associated to the kinetic energy and 
provides for the interaction between the magnetic field 
and the particles described by ip. Our model has 2 fixed 
parameters: b and the average density of particles, which 
is defined by 



1 



\ip\dS , 



(2) 



where J dS denotes the integral over the cross section of 
the cylindrical sample; this cross section will be taken as 
a disk and its radius as the unit of length. (If the radius 
is denoted by R, ip and A; scale as -R -1 , b and p as R~ 2 , 
and p as R°.) 

Our problem is to find the fields ip and Aj which mini- 
mize the thermodynamic potential J gdS for given b and 
p. For this purpose we minimize J(g — p,\ip\ 2 )dS, where 
p is a Lagrange multiplier. Variation of ip gives 



(iV - Ai - brO) 2 ip = pip 



(3) 



which is the usual time independent equation for a 
charged particle in a magnetic field. Here, instead of 
a single particle we consider a "condensate" ; by this we 
mean that instead of a single particle there are p par- 
ticles per unit volume, but they don't interact among 



themselves and they are all in the same state. The con- 
densate obeys the same equation as a single particle, but 
its current and magnetic influence will be proportional to 
p. In the limit p — ► 0, the thermodynamic potential ap- 
proaches irpp and Eq. (0) reduces to the Landau problem 



1- 



Variation of Aj gives 

VxVxAj = Re[tjj{iV - Ai - br§)ip] 



(4) 



where the bar denotes complex conjugation. This is just 
Ampere's law. If Ai (resp. ip) is fixed, then Eq. (0) (resp. 
([|)) is linear, but the system of both equations is nonlin- 
ear due to their mutual interaction. We choose a gauge 
such that Ai is parallel to 8. The boundary condition 
for Ai is continuity of V x Ai. For ip we will take the 
natural condition that its normal derivative vanishes at 
the boundary. 

There are solutions of the system (||)-(|1) with the axial 
symmetry of the problem. These have the form 



ip = K m (r)e~ 
Ai = A m (r)0 



(5) 



where to is an integer. For this form, (@)-(H) reduces to 
a system of ordinary differential equations. The value of 
to has to be chosen such that the minimum value of the 
thermodynamic potential is obtained. This value of to is 
an increasing function of b. Due to continuity, ip has to 
vanish along the axis of the sample for to ^ 0. 

We now ask whether there exist situations such that 
the minimizer of J gdS is not in the family (H). In the 
following, we will no longer consider p and b as inde- 
pendent parameters, but will focus on the value of the 
magnetic field for which the lowest value of J gdS among 
the solutions in the family (@) is shared by the wind- 
ing numbers to = and to = 1. For p < 10, it 
is found numerically that this approximately occurs at 
b = 1.924 + 0.171/9 + 0.0010V - 0.000036p 3 . We solve 
the problem in two stages: in the first stage we perform a 
variation in which ip has the form p + qr(l — r/2)e~ and 
A\e = a (r - 2r 2 /3) + ai(r 2 - 3r 3 /4) cos6». This varia- 
tion can be performed analytically, enabling us to obtain 
a qualitative picture of the minima and saddle points of 
JgdS. After we know what to look for, the system (0)- 
(||) is solved numerically. 

Besides the symmetric solutions with to = or m = 1 , 
we always find a nonsymmetric solution, but this has 
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a higher value of J gdS. (By "a solution" 
class of equivalent solutions.) However, for p > 2.7, a 
bifurcation occurs from the symmetric solution with m = 
1. This bifurcation is characterized by a migration of the 
nodal line away from the axis of the sample (Fig. nf). 
This behavior reminds of a transition numerically found 
for mesoscopic superconducting disks 01, except that in 
the present situation there is a single nodal line. 



tions in a broad range, and none of the trials lead to a 
lower potential. We might have missed a non-symmetric 
solution with lower potential, but this possibility would 
only strengthen our conclusion that axial symmetry is 
broken. 
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FIG. 1. Distance of the nodal line, r no do, from the axis of 
the sample (divided by the sample radius). For a solution in 
the family (pi) with m = 1, the nodal line is at the axis of the 
sample. At p ~ 2.7, a new (asymmetric) solution bifurcates 
from it. 

We find that this nonsymmetric solution has lower 
thermodynamical potential than any couple of fields in 
the family (pft (Fig. 0). This means that the minimizer 
of J gdS is not in the symmetric family rtq) , thus giving 
rise to symmetry breaking. 




FIG. 2. Thermodynamic potential of asymmetric solutions, 
in comparison to those of symmetric solutions. <? S ym is the po- 
tential density that gives the lowest integral J g Byln dS in the 
symmetric family with m = or m = 1. (For the magnetic 
field considered, m — and m = 1 give the same integral.) 
For p > 2.7, the lowest potential is obtained for an asymmet- 
ric solution. 



One might suspect that the solutions of @-(Q) we have 
found are not those with the lowest thermodynamical po- 
tential. This is unlikely. For p — > 0, our symmetric so- 
lutions approach well known analytic results; for p = 10 
we have solved (||)-(||) for initial symmetric configura- 



